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Abstract 



We perform a comparison between different lattice regularizations 
of the Dirac operator for massless fermions in the framework of the sin- 
gle and two flavor Schwinger model. We consider a) the Wilson-Dirac 
. operator at the critical value of the hopping parameter; b) Neuberger's 

overlap operator; c) the fixed point operator. We test chiral properties 
of the spectrum, dispersion relations and rotational invariance of the 
mesonic bound state propagators. 
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1 Motivation and Introduction 



The Nielsen-Ninomiya theorem and the Ginsparg- Wilson j2] condition 
(GWC) provide us with the crucial information under which circumstances 
[3j remnants of chiral symmetry may stay with a lattice action for fermions, 
like it is e.g. the case for overlap fermions jl] or fixed point actions [5]. 
Recently Liischer [S| has pointed out the explicit form of the underlying 
symmetry and indicated possible generalizations. 

Fermion actions for massless fermions satisfying the GWC 

i{l?, 7 5 } = a (V^RV) , (1) 

where T> is the lattice Dirac operator, violate chiral symmetry up to a local 
term 0(a); the r.h.s. of the above equation is local since R is. 

It was pointed out in [3| that actions which are fixed points under real 
space renormalization group (block spin) transformations (BST), are solu- 
tions of the GWC. R is then local and bounded and as a consequence the 
spectrum of T> in complex space is confined between two circles 

I A Train] — ^min > |^ T'rnaxl — ^max i (2) 

where the real numbers r min and r max are related to the maximum and mini- 
mum eigenvalue of R respectively. For non-overlapping BSTs R = \ and (J2J) 
reduces to | A — 1 1 = 1, i.e. the spectrum lies on circle. 

Independent solutions of the GWC are provided by the overlap formalism 
0], which allows the formulation of chiral fermions on the lattice. These 
solutions are obtained in an elegant way, as shown recently by Neuberger 
[II IE], through some projection of the Wilson operator with negative fermion 
mass. Even in this case we have R = \ and |A — 1| = 1. 

For the Schwinger model (2D QED) we are in a situation where we have 
access to three different lattice Dirac operators for massless fermions, namely 
the original Wilson operator D Wi at k c (/9), the Neuberger-projected operator 
Pnc and a numerically determined and therefore approximate fixed point 
operator £>f p HDj- Studying these alternatives we may ask the following 
questions: 

• For given gauge field configurations: What is the relation between real 
eigenvalues of the lattice Dirac operator and the geometric topological 
charge? To what extent is the Atiyah-Singer Index Theorem (ASIT) 
realized in these lattice environments? 
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• In the continuum, the eigenvalue distribution is related to the conden- 
sate (ipip) (Banks-Casher formula JI]). What about the lattice theory 
for the given Dirac operators? 

• The first two questions concern chiral symmetry and the phenomenon 
of fermion condensation. Important for the eventual study of the con- 
tinuum limit of the full theory are also spectral properties: What is the 
behavior of e.g. dispersion relations? Concerning off-shell properties, 
what about recovery of rotational invariance of the propagators? 

Here we perform a Monte Carlo simulation for the one- and two-flavor 
Schwinger model with gauge fields in the compact representation and the 
three different lattice Dirac operators. For the two-flavor model one expects 
a massless bound state in the chiral limit, which is of particular interest in 
this framework. We stress that the ensemble of gauge configurations used 
is the same in all three cases, the sampling being performed using the one- 
plaquette standard gauge action. The unquenching is obtained through the 
multiplication by the fermion determinant. 

2 Lattice Dirac Operators 

The three fermion actions may be written tpT>tp, where the fields at each site 
are two-components Grassmann variables ■?/>, ip and T> is a matrix in Euclidean 
and Dirac space (lattice Dirac operator). For two flavors the number of fields 
duplicate and the action becomes u T> u+dV d (with independent Grassmann 
fields u, u, d, d). 

All three Dirac operators are non-hermitian but have 75-hermiticity: 
75D75 = T>>. Their eigenvalue spectrum is therefore symmetric with regard 
to complex conjugation. 

2.1 Wilson Dirac operator 

We write the Wilson Dirac operator in the form 

£> W i (x, y) = (m + D) l xy - - ^ [(1 + °>) U xy + (1 - c M ) U\ x S x>y+fM . 

(3) 

The hopping parameter k is related to the (bare) quark mass m through the 
relation k = 1/ (2m + 2D). 
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The chiral limit is obtained in this environment for k — > K c (/3), where 
k c (/3) should be determined in some way (see the following). Thus we will 
work at D — 2, and k = k c {(3) (corresponding massless quarks) as discussed 
below. 

2.2 Neuberger's operator 

Neuberger suggests [7j to start with the Wilson Dirac operator at some value 
of m G (-1,0) corresponding to ^ < K < 2£) 1 2 and then construct 

D Nc = l + 75 e(7 5 D Wi ) . (4) 

We call the actual value of k used in the above definition K Ne . Some words 
about the choice of K Ne : according to [8 it is arbitrary, in the sense that any 
(strictly negative) value of m in the interval (—1,0) reproduces the correct 
continuum theory (see also the discussion on a suitable choice of Kn e in [T2]). 
but it may be optimized with regard to its scale dependence by looking for 
example at the behavior of the (projected) spectrum. Comparing expectation 
values of operators like (ipip) for different /% e one has to take care of the 
proper normalization ^3] • One may see from the comparison with free lattice 
fermions that there is a (trivial) factor of y/m, i.e. (^■0) = (V^)Ne/|^| in 
our conventional Whenever not mentioned otherwise we choose K^e = \ 
(m = —1) for our exploration; we also discuss results with some smaller 
values. 

The operative definition of e( 7 5 £>wi ) entering the above equation is: 

e ( 75 V Wl ) = U Sign(A) U ] with 75 V Wl = U A U ] . (5) 

Here Sign(A) denotes the diagonal matrix containing the signs of the eigen- 
value matrix A obtained through the unitary transformation U of the hermi- 
tian matrix 75 Pwi • There are various efficient ways to numerically find 
£>Ne without passing through the diagonalization problem, prohibitive for 
D = 4 [E] (for D = 2, see also J3]). In our simple context computer time 
is no real obstacle and therefore we use the direct definition ©, explicitly 
performing the diagonalization. 

As observed in jH] , gauge configurations with non-zero topological charge 
imply exact zero eigenvalues of V Ne ■ The subsequent inversion, necessary to 

1 We thank H. Neuberger for pointing this out to us. 
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find the quark propagators, is then not possible. The correct way to proceed 
is to introduce a regulator cutoff-mass: D — > T> + /x 1 and then consider the 
limit jj, — > 0. 

2.3 Fixed point Dirac operator 

In |H| the fixed point Dirac operator was parameterized as 

3 

Thv(x,y) = Pi(f) a i u i x J) > withy = x + 5f. (6) 

Here / denotes a closed loop through x or a path from the lattice site x to 
y = x + 5f (distance vector Sf) and U (x, f) is the parallel transporter along 
this path. The <Tj-matrices denote the Pauli matrices for % = 1,2,3 and the 
unit matrix for 2 = 0. The action obeys the usual symmetries as discussed 
in altogether it has 429 terms per site. The action was determined for 
gauge fields distributed according to the non-compact formulation with the 
Gaussian measure. There excellent scaling properties, rotational invariance 
and continuum-like dispersion relations were observed at various values of 
the gauge coupling (3. 

In ^U] the action was studied both, for compact and the original non- 
compact gauge field distributions. In the compact case the action is not 
expected to exactly reproduce the fixed point of the corresponding BST, but 
nevertheless it is still a solution of the GWC; violations of the GWC are 
instead introduced by the parameterization procedure, which cuts off the 
less local couplings. We demonstrated that indeed the spectrum is close to 
circular, somewhat fuzzy at small values of /3 < 2 but excellently living up 
to the theoretical expectations of jHj at large gauge couplings (3 > 4. 

Here we study the action only for the compact gauge field distributions 
in order to allow a direct comparison with the other lattice Dirac operators. 

3 Simulation Details 

Uncorrelated gauge configurations have been generated in the quenched setup. 
However, we are including the fermionic determinant in the observables: 
all the results presented here are obtained with the correct determinant 
(squared, for two flavors) weight. From earlier experience [HI E3] we know 
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that this is justifiable for the presented statistics. We perform our investiga- 
tion on three sets of 5000-10000 configurations at /3 — 2, 4 and 6. 
The so-called "geometric definition" of the topological charge is 



we keep track of its value for all our configurations. The configurations have 
been well separated by 3r int , the autocorrelation length for Q G . 

For P Wi we need to determine k c ((3). We use PCAC techniques for this 
purpose JH] determining k c for the unquenched 2-flavor case. 

For each configuration we then build XVi (at k c ((3)), V-^ e a.nd Vp p as dis- 
cussed. Each lattice Dirac matrix is diagonalized to obtain the complex 
eigenvalue spectrum. This is somewhat time-consuming due to the non- 
hermiticity. Furthermore the inverse (the quark propagator) is determined. 

In the 2-flavor Schwinger model one expects J7j (for a recent discussion 
cf. PHI) one massive mode (called rj by analogy) and a massless flavor-triplet 
(called 7i"). The corresponding momentum-projected operators are 



r](p,t) = ^ (u(x l7 t) (Jxuix^t) + d(x 1 ,t)a 1 d(x 1 ,tf) , (8) 



K3(jP,t) = e%PXl ( M ( x i> t) °~i u(xi, t) - d(xi, t) <7\ d(xi,t)j . (9) 



Their correlation functions define by their exponential decay the correspond- 
ing energy functions E(p) and thereby the dispersion relation. In the 1-flavor 
case only the massive mode is there. 

We also study rotational symmetry via the correlation function 



measured for all 2-point separations. 

As mentioned before, in order to avoid numerical problems with inversions 
for the Dirac operator with (almost or exact) zero eigenvalues we introduce 
a small regulator mass //. It turned out that for \i = 0(1O~ 3 ) or smaller the 
result is practically insensitive to this cut-off, the inversion algorithm still 
working properly. 






(10) 
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4 Discussion of the Results 



From the definition (jlj) of V Ne one could naively expect that its spectrum is 
obtained from the starting Wilson operator at a proper value of K^e by simply 
projecting the eigenvalues onto the circle on the complex plane |A — 1| = 1. 
This is not quite the case, although it becomes more and more so for larger (3 
approaching the continuum limit. In particular the real modes are projected 
either onto A = or to A = 2 (see fig^). In the case of a pair of real modes of 
opposite chirality, they split into two conjugated complex eigenvalues on the 
circle. Finally one is left with zero modes of only one definite chirality - in 
agreement with the so-called vanishing theorem valid in D=2 ^Hj _ and an 
equal number of modes A = 2 with the opposite chirality. As a consequence 
of the process of splitting of chirality pairs, the number of zero modes in 
£>Ne is smaller than those of real small eigenvalues of £Vi • 




-1.0 0.0 1.0 -1.0 0.0 1.0 

Re Re 

Figure 1: Eigenvalues of V Wi (open circles) and of the term 756(75!^) 
(full diamonds) in V^ e , left: KNe = §, right: K^e = §• 

In fig. [U we compare - for a fixed gauge configuration (at (3 — 6) - the 
eigenvalues of £>wi (at « = 5 and ~; n.b. this is still above k c ) with those 
of V Nc — 1 = 756(75^^^1) resulting from the projection (jlj). We find that 
the negative real eigenvalues of Dwiare projected to —1 (corresponding to 
A = 0): their number, counted according to the signs of their chirality, agrees 
with the number of zero modes of . In particular gauge configurations 
for smaller (3 have more eigenvalues on the real axis, with a distribution 
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Dirac op. 


p((3 = 2) 


p(P = 4) 


p(P = 6) 


£>Wi 


74.22 


99.70 


100.00 




74.22 


99.70 


100.00 


V Fp 


96.58 


100.00 


100.00 


£>wi 


62.56 


97.36 


99.78 




100.00 


100.00 


100.00 




91.00 


99.84 


99.96 



Table 1: Upper part: percentage p(/3) of configurations where the Atiyah- 
Singer index theorem is fulfilled. Lower part: percentage p(/3) of configura- 
tions where the vanishing theorem is fulfilled. 

density becoming broader when decreasing (3. Increasing (decreasing) 
the spectrum of £\yi would shift towards left (right), and so more (less) 
zero modes might be obtained as a result of the projection. For small j3 
there is no clear distinction for £\yi between the physical branch of the real 
spectrum and the eigenvalues due to doubler modes. This uncertainty is 
discussed in the framework of the overlap method [20] and the Wilson- and 
Sheikholeslami-Wohlert action |21j . 

For £> Ne , as for D Wi , not all zero modes can be equivalenced to the 
geometrically defined topological charge of the gauge configuration. Of course 
one can still define the topological charge as the number of zero modes, or 
equivalently the number of A = 2 modes, counting properly the chirality; 
in the case of the fixed point action one obtains in this way the fixed point 
topological charge of the lattice configuration 5J. This ambiguity vanishes 
towards larger /3, approaching the continuum limit. Table ^ (upper part) 
shows the percentage p(/3) of configurations where the number of zero modes 
counted according to the sign of their chirality (jir — n£) (cf. the discussion 
in [22]) a g rees with the geometric topological charge. The agreement of the 
first two lines is trivially explained, since the real modes counted for X\yi are 
those projected to A = (or split into complex pairs in the case of a real pair 
of opposite chirality) for T>^ e . This agreement stems from the choice of /% e 
as "cut-off" for the counting of the real modes of £\yi • 

In tabled (lower part) we confirm that for DNethe zero modes have just 
one definite chirality, whereas for the other actions both chiralities contribute, 
displaying a violation of the vanishing theorem (in the case of the V Fp we 
believe that this violation is an effect of the truncation of less local couplings). 
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For larger f3 all three actions recover the vanishing theorem and the ASIT: 
(riR - n L ) = Q G . 

The density distribution of eigenvalues on (for T>-^ e ) or almost on (for 
Ppp ) the circle agrees with each other at small eigenvalues, with improving 
agreement for increasing (3. For « Ne = | already at j3 — 6 the densities are 
indistinguishable within the statistical errors for Im(A) < 0.6 (Fig. here 
we adopt the definition in (TQ3 f° r the projection of the density distribution 
of eigenvalues to the imaginary axis). Choosing another will produce 
another eigenvalue distribution. Comparing these with each other and with 
the value of ^1] one has to take care of the proper normalization of the 
fermion fields, as discussed above in sec. 12.21 



0.150 
0.125 
0.100 
0.075 
0.050 
0.025 
0.000 



-3 



-2 



-1 



Figure 2: The unquenched eigenvalue density distribution projected from 
the circle onto the imaginary axis at /3 = 6 for T>p p (thick lines) and X>nc 
with «Ne — \ (thin lines). The horizontal line denotes the continuum value 
at infinite volume. 

A suitable lattice definition of (ipip) was suggested in the framework of 
the overlap action [22] , which can be generalized for any fermions obeyed the 
GWC [21] (see [TU] for an application). Also this quantity shows excellent 
agreement for at = \ and Pfp • We obtain for the two actions (lattice 
units): 0.073(4) and 0.072(7) respectively for (3 = 4, 0.062(3) and 0.063(3) 
for /3 — 6 (to be compared to the continuum values 0.080 and 0.065 |25j). 
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Choosing = | with the appropriate normalization (with the factor for 
free fermions) we find 0.058(1) at f3 — 6. The condensate has been studied in 
the 1-flavor Schwinger model in the overlap formalism already in [20] , where 
also good agreement with the expected values has been demonstrated. 




1.5 
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0.5 



0.0 
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Figure 3: The dispersion relation E(p) for (upper figure) the ir and (lower 
figure) the rj propagator for the three actions studied (lattice size 16 2 , (3 = 6); 
squares: Wilson action, diamonds: Neuberger action, circles: fixed point 
action. 

Discussing the mass spectrum, we find that for V^ e (n^ie = 5) ^he van " 
ishing of the 7r-mass is excellently realized (see fig. this is an amazing im- 
provement with respect to "D-yvi > where the tuning to k c creates non-irrelevant 
technical problems. The behavior of the energy dispersion for non-zero mo- 
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menta shows no improvement compared to £\yi > while for Vp p it is almost 
linear as expected for a massless particle in the continuum: Vp p eliminates 
the cut-off effects producing a continuum- like propagator for all momenta. 
The massive state displays qualitatively similar behavior. 

Decreasing K Ne in the definition of V Ne from the suggested value | down 
to | somewhat enlarged the values of E(p) slightly above those of XV; , but 
still similar in the overall shape. Changing further to even smaller K^e closer 
to \ we observe, that the propagators do not reach asymptotic behavior on 
the lattice size studied (even for K^ e > k c ) and no mass plateaus can be 
identified. We suspect, that this indicates larger corrections to scaling. 

We also studied the correlation function for the three actions. The 
rotational symmetry properties for T>-^ e are comparable to those of Pwi • Ac- 
tion Dpp shows the best rotational invariance. In summary, we find that the 
real-space correlations functions and the dispersion relations are not notice- 
ably improved by T>^ e , as compared to those for Vp p , which show a behavior 
substantially closer to the continuum. 

It is expected that £>N e is automatically 0(a) corrected |2Zj and thus 
improves scaling for the on-shell quantities; without at least introducing im- 
provement of the current operators one would not expect improvement for 
the propagators as exhibited by our results. 

We conclude that the identification of zero modes with the geometric 
topological charge for T>^ e agrees with that for XVi > if the real modes are 
counted according to their chirality. T>p p shows generally better behavior. 
The vanishing theorem, however, is satisfied automatically for T>^ e , since 
only zero modes of one chirality result from the projection. The bound state 
masses and condensate come out similarly. The rotational invariance and the 
dispersion relations of the current propagators are not significantly improved 
for 2?Ne as compared to £>wi , but they are definitely better for Vp p . 
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